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Abstract 

Non-Abelian flux-tube (string) solutions carrying global currents are found in the 
bosonic sector of four-dimensional J\f = 2 super-symmetric gauge theories. The specific 
model considered here possesses U(2)i oca ixSU(2) g i 0 b a i symmetry, with two scalar dou¬ 
blets in the fundamental representation of SU(2). We construct string solutions that are 
stationary and translationally symmetric along the x 3 direction, and they are character¬ 
ized by a matrix phase between the two doublets, referred to as “twist”. Consequently, 
twisted strings have nonzero (global) charge, momentum, and in some cases even angular 
momentum per unit length. The planar cross section of a twisted string corresponds to a 
rotationally symmetric, charged non-Abelian vortex, satisfying first order Bogomolny-type 
equations and second order Gauss-constraints. Interestingly, depending on the nature of 
the matrix phase, some of these solutions even break cylindrical symmetry in M 3 . Al¬ 
though twisted vortices have higher energy than the untwisted ones, they are expected to 
be linearly stable since one can maintain their charge (or twist) fixed with respect to small 
perturbations. 


1 Introduction 

Vortex and string-type solutions appear in many models, and as they have many applications 
there is an enduring interest for them. In the plane, a vortex corresponds to a cross-section 
of a straight string in three spatial dimensions in the plane orthogonal to its direction. In 
the context of spontaneously broken gauge field theories by scalar fields, the paradigm is the 
Abrikosov-Nielsen-Olesen (ANO) vortex [I] associated to the breaking of an Abelian gauge 
group. The ANO vortex corresponds to the planar cross-section of a static straight, infinitely 
long magnetic flux-tube, with quantized magnetic flux and SO(2) cylindrical symmetry . ANO 
vortices have an (integer) winding number, proportional to their quantized magnetic flux, which 
is also responsible for their stability. For a fixed winding number, ANO solutions form a one- 
parameter family, depending on the mass ratio f3 = m s /m v , where m s , resp. m v denote the 
mass of the scalar resp. of the vector field. In the special case f3 = 1, the energy of a vortex is 
proportional to its winding number and surprisingly vortices of like fluxes do not interact 
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[2] , For this special value of the coupling, minimal energy vortices satisfy a set of first order - 
Bogomolny - equations, which are easier to solve than the held equations m la¬ 
in non-Abelian gauge theories with or without a Chern-Simons term, non-Abelian vortices 

(some of them with an electric charge) were first obtained in Ref. (3|. Since the seminal papers 

[3] started the investigation of vortex-string solutions in supersymmetric non-Abelian gauge 
theories, the subject continues to attract attention. A simple model containing the essential 
features is an U(N C ) gauge theory, coupled to Nf scalar fields in the fundamental representation, 
where Nf > N c . Vortex solutions in such theories are usually referred to as non-Abelian vortices, 
(NAVs). NAVs have attracted considerable interest, since they are at the heart of intriguing re¬ 
lationships between two-dimensional sigma-models and four-dimensional gauge theories. NAVs 
in a U(N) gauge theory possess “orientational moduli” whose low-energy dynamics is described 
by versions of CP^ -1 sigma-models on the string worldsheet. Moreover, NAVs also held the 
promise to be relevant to bring us closer to a description of quark confinement; for a review, see 
[6] and references therein. Static NAVs are absolute minima of the energy functional in a fixed 
topological sector, characterized by the winding numbers associated to the Cartan subalgebra 
of the gauge group. Moreover NAVs satisfy first order, Bogomolny-type equations admitting 
rather complicated, static multi-vortex solutions. Most remarkably a complete description of 
the NAV moduli space has been found m m- An interesting application of non-Abelian vor¬ 
tices is illustrated by confined monopoles emerging as junctions of NAVs with different moduli 
P HU]- The dynamics of NAVs based on the moduli approximation has been worked out in 
Refs, fm, K, H3]. In addition to NAVs and related monopoles, domain walls have also received 
due attention nail, and a number of intriguing relations between moduli spaces of monopoles 
and domain walls have been discovered |15j . 

In this paper we point out that by allowing for space-time-dependent phases among the 
scalar fields, which we shall call “twisting”, new families of charged vortex-strings arise. We shall 
restrict our attention to a U(2) gauge group broken by two scalar doublets with an appropriate 
scalar potential compatible with super-symmetry (SUSY), but a generalization to other groups 
should not be too difficult. It is known, that imposing the usual space-time symmetries on 
held configurations leading to vortex-type solutions, such as translational invariance in time 
and along the, say, x 3 or z direction, scalar fields may have a phase, with a linear dependence 
on (t,z) [[[SI CEEJ- In general, the energy of stationary configurations with a nontrivial 
(t , z)-dependent phase is bigger than that of the static ones, however this does not make them 
necessarily unstable [16]. 

In this paper we systematically investigate straight vortex-string solutions in the simplest 
theory admitting non-Abelian vortices, when the scalar fields, $, possess a (t , z)-dependent 
phase, i.e. 

$(/) = (I>(x r )e^ MuJaXa , i = 1, 2, a = 0, 3 , 

where the flavours correspond to the columns of the (2x2) matrix $, M is a constant Hermitian 
matrix acting on the flavour indices, and uj a is vector in (t , z) plane. As it is well known 
da ei fl9] when tu a is a light-like vector , the held equations in the (a; 1 , x 2 ) plane decouple 
completely from those in the (t , z) directions. The chromo-electric components in the (t , z) 
plane are determined by a set of linear, second order equations, (Gauss constraints), depending 
on the solution in the (x 1 , x 2 ) plane. Simple analysis of the Gauss constraints makes it very 
plausible that any NAV can be twisted by an arbitrary twisting matrix, although we have not 
attempted to formally prove this. There are global currents flowing in the z direction of an M- 
twisted vortex-string which has nonzero charge, momentum and, unless M is specially aligned 
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in internal space, angular momentum. A twisted NAV has some genuine three-dimensional 
structure, as there are global currents flowing along its axis, the fields have some nontrivial 
components orthogonal to its symmetry plane, hence it is more conveniently thought of as a 
twisted string. 

The simplest twisted NAV can be obtained from the rotationally symmetric “elementary” 
vortex solution of Ref. |5], in which case, for a general twisting matrix the problem reduces 
to a single second-order Gauss-constraint, which can be easily analyzed. In fact our twisted 
NAV solution turns out to be the same object as the “dyonic” vortex solution found in a 
”mass-deformed” SUSY gauge-field theory G2t The bulk of our paper concerns the twisting of 
“composite-coincident” NAVs p], corresponding to the rotationally symmetric superimposition 
of elementary vortices. We characterize these superimposed NAVs by a relative winding number 
and a moduli parameter. Interestingly we find that when the twisting matrix M, contains off- 
diagonal components the twisted string looses cylindrical symmetry in three dimensions. 

The energy of an M-twisted string is higher than that of an untwisted one, it is given as 
a sum of the usual “magnetic” energy per unit length proportional to the magnetic flux, and 
of an “electric” contribution due to the rotating phase. The actual magnitude of the “electric” 
contribution depends essentially on the components of the twisting matrix, M. When M is 
diagonal the magnitude of the electric energy of a twisted NAV is typically much smaller than 
its magnetic one, in fact, it can be made arbitrarily small when the untwisted NAV is sufficiently 
close to a diagonal one. In the case when M contains off-diagonal components, the magnitude 
of the electric energy is comparable to the magnetic one. 

It seems to us that all the above-mentioned properties make twisted NA vortex-strings of 
some interest and worthy of further investigations. 

The plan of the paper is the following. We introduce in Sec. [2] the theory that we shall study 
and proceed to the dimensional reduction by splitting four-dimensional Minkowski space into 
planar and temporal-longitudinal coordinates. Making an appropriate Ansatz we get a dimen¬ 
sionally reduced Lagrangian. Assuming that coupling constants satisfy Bogomolny conditions 
we present the minimal energy first-order equations satisfied by untwisted solutions and write 
the energy as a sum of electric and magnetic contributions. We discuss how the rotational 
symmetry of the Ansatz, or the absence thereof, depends on the properties of the twisting. In 
Sec. [3]we consider twisted vortices, both elementary (Sec. 13.11) and composite (Sec. 13.21) ones, 
presenting the numerical study of the solutions. In Section [4] we summarize the results and 
present our conclusions. Finally our notations and conventions are given in Appendix [Aj and 
some numerical data is given in Appendix iBl 

2 Dimensional reduction 

The 3+1-dimensional theory we consider is defined by the following Lagrangian: 

C = ~F^ - + Tr (D^D^ -V, (1) 

^9 i ^9 2 

where is the Abelian field strength tensor and G" tv (a = 1, 2, 3) is the non-Abelian, SU( 2) 
one. The two scalar doublets are encoded in the matrix (‘Fjm with i = 1,2 being the gauge 
(“colour”) and A = 1,2 the “flavour” index; the trace is taken over the flavour indices; the scalar 
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potential, V — V\ + V 2 , can be written more explicitly as 


( 2 ) 

(3) 


Vi = 


v 2 = 


y(Tr$ f $ - 20 2 = y (& a $a ~ 2 ^ 
y (Tr $V “$) 2 = ^{*%< 7%* jA ) 2 = y 


($l$i - 


$^$ 2 ) 2 + 4|^$ 1 



where summation over repeated indices is understood, except when the contrary is indicated. 
For more details on notations and conventions, see Appendix [A] 

The fields transform under the 17(2) gauge symmetry as 


<D^exp(*A)17<h, C^UC^ + 2iUd^, -> A^ + 2d^A , (4) 


where U(x) G SU( 2), A(x) is a real function. The flavor symmetry acts from the right on the 
scalars as 

V G SU{2). (5) 

Let us now consider stationary and translationally symmetric fields in the x 3 direction. It 
will turn out to be convenient to split four-dimensional Minkowski coordinates as x M = (x a , x l ) 
with a — 0, 3, % — 1, 2. Since the symmetries are generated by two commuting vector fields, 
there exists a gauge where the symmetric gauge fields are simply independent of the coordinates 
x a [20] . At the same time, however, the scalar doublets, <3>, may still depend linearly on x a 
through an SU(2) phase, i.e. the most general symmetric Ansatz can be written as 


A^x v ) = (A^x- 7 ), A a (x J )), (6a) 

C;(xn = (C?(x>) : C a a (xi)), (6b) 

$(x 1 ') = $(x i )e3 Mw “ iea , (6c) 

where M is a constant Hermitian matrix, 

M = m h a a = m° 1 + mV a , (7) 


c o a is a vector in the (t,z) plane. Straight flux-tube/string solutions of the theory defined 
by Eq. ([T]), described by the Ansatz (J6]) with a non-trivial M, will be referred to as ‘twisted’ 
strings. 

The form of the Ansatz, Eq. ([6]) also restricts the symmetries of the model. Those sym¬ 
metries, which preserve the Ansatz are flavour transformations which commute with the twist 
matrix M, 

$ -> $17 , V = exp (iMS) , (8) 

and all gauge transformations where Ud a U' and d a A only depend on x l and not on x“. The 
flavor-current generating the transformations (J8]) can be written as 

K, = m h K J , (9) 

with K* defined in Eq. (l57l) . 

Being Hermitian, the twist matrix, M, can always be diagonalized by a unitary matrix, Vm, 
such that M = VmMdVm\ with Mo being diagonal, i.e. 

^(x i )ei MLJaX °‘ = ^)V M e^ MDUaafX V M ^ ■ ( 10 ) 
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Therefore, twisting a configuration with the matrix M, is equivalent to twisting a suitably 
flavour-transformed (with Vm) configuration with the diagonal matrix, M D . As the theory 
considered is SU(2) flavour symmetric it makes no difference which form of twisting one uses. 
In what follows we shall work with the non-transformed twisting matrix, M in Eq. (Iflcl) . 

The Ansatz (jHD yields the dimensionally reduced Lagrangian 

£ + ^i d - A ° a ' A ° ~ 4ii< G “T + k D ' clDiC “ 

- jklfCJC"”) 2 - (C‘Cp)(C k ‘C bl> )] (11) 

^9 2 

- Tr(A$) f + 2 Tr - C a $) t (A$M - C“4>) - V , 

where C a = A a + C£a a . It is now convenient to introduce basis vectors in the x a plane, (uj a , u: a ) 
such that u> a Cu a = 0 if oj 2 = LO a u> a ^ 0, while for the light-like case ( u 2 = Co 2 = 0) uj oi A ol ^ 0. We 
remark that in Eq. (fTTl) only the a indices are raised or lowered by the induced Minkowskian 
metric in the (x°,x 3 ) plane, and the repeated i,j -type lower indices are summed with (+,+) 
signature. It turns out that the dimensionally reduced Lagrangian (fill) , is closely related to the 
trivial reduction from four to two dimensions of the “mass deformed” theory considered in Ref. 
[12] . In the absence of twist, M = 0, (fTTl) corresponds to the trivial reduction of the theory 
to two dimensions, whose solutions are the non-Abelian vortices discussed in detail in Ref. [8]. 
With respect to the c o ,u> basis the gauge field components in the (x°,x 3 ) plane are expressed 
as 

A a = Ld a A + id a A , Cq = Ld a O a + Ld a C a . (12) 

The field equations can be grouped according to variations with respect to A , C resp. A , C. 
Consider first the variational equations with respect to A,C: 

AA — y Tr($ f C$), 

A: A ° a = -f ®A(Ao a + C a )$ A ~ A a , 

uj 2 C b (C b C a - C b C a ) if (U V 0 
2(uu) C b {C b C a - C b C a ) if a; 2 = w 2 = 0 , 

where A = didi with Euclidean metric summation. From Eqs. (TT3I) we obtain the following 
integral identity: 



(13a) 

(13b) 

(13c) 


d 2 x 


- -4 aA) 2 - 4(^A a ) 2 - I^a^a + \c a A° 

2 9i g z 2 2 c /2 


= 0, 


(14) 


where A = A (A 2 /gf + C a C a /g 2 ). Assuming finite-energy boundary conditions and global 
regularity, the integral of A is zero. For the case when u is light-like, oo 2 = 0, one finds 
C a A a = 0, therefore Eq. (TH1) enforces A = 0, C = 0. When cu 2 < 0, Eq. (fTTl) implies once 
more the vanishing of A and C, since then C a A a < 0. In the case of a space-like oj vector, 
oj 2 > 0, Eq. (TT4l) is not sufficient to exclude the existence of non-trivial solutions of the Gauss 
constraints, (fT3j) . It is consistent, however, to assume A = 0, C = 0, even for u> 2 > 0, since A, 
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C satisfy homogenous equations. Assuming A = 0, C = 0 the remaining held equations are 

(15a) 


= - Tr [($ f C - , 

9i 2 


-4= - Tr + (A - M)$V) $] , 

92 ’ 


1 


2 

i 




1 DjG” = - Tr |>[> t rT' ; /),'[> - /i/l>V‘.[, - ‘^ e °-t*C b D J C ! 


92 


D ‘ D ‘* = »• 


92 

8V lo 2 

- — [($M - 2C<h)M + C 2 $] 


(15b) 

(15c) 

(15d) 

(15e) 


The total energy of an M-twisted string can be written as the sum of an “electric” and of a 
“magnetic” part: 

E = I d 2 xT° 0 = I{S 0 + A) d 2 x =Eo + E x . (16) 

where the “electric”, resp. “magnetic” densities, S 0 , S\ are defined as 


So = GaiGli + TT^EaiFai + (D a $ A )'D a $ A 
^9 i 


2 9l 


— (c^o + col) 


1 ( diA) 2 + 9-{DiC a ) 2 + \ Tr($Af - C<h) t ($M - C$) 


L2^f 


2<? 2 2 


1 


A = tt (^) + tt(^) + IA<M 2 + v. 
^#2 ^1 


(17) 


(18) 


A straightforward computation shows that using the Gauss constraints, Eqs. (115af) - f115bjl . the 
“electric” density, S 0 can be expressed as 

1 


^0 — "j^O + W 3)(A + Q) , 

where A = A {A 2 /g\ + C a C a Ig\ ) and 

Q = Tr [$t($Af - C$)Af] . 


(19) 


( 20 ) 


Note that the (t, z)-components of the flavor current in Eq. (JUJ) can also be expressed in terms 
of co and Q as 

K a = -~^Q ■ (21) 

In most work on non-Abelian vortices the supersymmetry-induced relations between the cou¬ 
plings, Ai = gf, A 2 = g%, have been assumed. These relations ensure that S\ can be expressed 
as a sum of squares and a topological term, leading to first-order, Bogomolny-type equations in 
the (x 1 , x 2 ) plane. Then minimal-energy, untwisted solutions of the second-order held equations 
(115cD (115eH are obtained by solving the following hrst-order equations: 


T.J.. = +A, t 'Tr'I't.I. _2f), 

Gi = : FfealtaVS, 

D& = Tiei k D k $ , 


(22a) 

(22b) 

(22c) 
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while for solutions of Eqs. (1221) the “magnetic” energy density simplifies to 

£\ = ±^e ik F ik T ie ik di • (23) 

Equation (1231) implies that the total “magnetic” energy, Ei, is given by the net Abelian flux 
through the (x^x 2 ) plane. Let us quote here the actual value of the “magnetic energy” in Eq. 
(l23lh Ei, for configurations considered in this paper, characterized by winding numbers tla , rriA 
as in Eq. (1271 and subject to the first-order Eqs. (|22l) : 


Ei = 27r£|n 1 + m 2 \ = 27r£|n 1 + mi + N \, 


(24) 


which of course also holds for the Ansatz (IHTTh In the case when oo a is light-like, the held 
equations in the (x^x 2 ) plane, Eqs. (j!5c[l (I15e[) . decouple from the Gauss-type constraints, 
Eqs. (I15ap - fll5b|) . and become identical to those corresponding to untwisted vortices. Therefore 
the problem of Ending twisted non-Abelian strings for co 2 = 0 reduces to solving Eqs. (j 15 all 
(U5bj) in the background of a non-Abelian vortex in the (x^x 2 ) plane. In the present paper 
we shall consider a light-like u vector, and concentrate on twisting vortex solutions of minimal 
energy satisfying the Erst-order Eqs. (1221) . It is left for future work to clarify if for to 2 ^ 0 there 
exists solutions analogous to the twisted vortices of Refs. mm- 

An M-twisted string has momentum flowing along the z direction, which is easily obtained 
from the stress-energy tensor. The longitudinal momentum, P = Tq 3 , carried by a twisted 
string, can be recast exploiting the Gauss constraints, Eqs. (II5aj) (115bh . as 


P — T 03 — -UqU3(A. + Q) 


(25) 


M-twisted strings may also have angular momentum, J, as it can be seen from the angular 
momentum density 

J — Tqd = —Fo r Ftf r -|—2 A Tr F)q& + Tr ; (26) 

9 1 9 2 

however, to compute J in a more explicit form, one needs a parametrization of the angle 
dependence of the fields. This will be presented in the next section. 


2.1 The Ansatz; rotational symmetry and its loss 

Let us now impose rotational symmetry in the x^x 2 plane to the fields. Denoting the usual 
polar coordinates in the plane as x 1 = rcos d, x 2 = r sin d, rotational symmetry implies that 
by 'd-dependent gauge transformations one can always achieve 


<94-V = 0 , $a(x 1 ) = {exp(inA’d)®Ai(r),exp(im A d)$A2(,r)) , A = 1,2. (27) 

In order to ensure consistency with the U (2) gauge and the global SU (2) flavour symmetry, the 
integers, uatIxia-, satisfy the following relation: n 2 — n\ = m 2 — mi = N, which can also be 
expressed on the two scalar doublets as 


<f>(x*) 


k(j)i (r)e* ni1? ijji(r)e mi ' & \ (\ 0 \ 

\cj) 2 (r)e im ^ 'i/j 2 (r)e imi ' & J 1^0 e im ) 


$o(x*)e 


iNi? 

5 


(28) 
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where N = Diag{0,lV}, which encodes the relative winding between the two flavours. Since 
in general the twisting matrix, M, mixes the two flavours, when N ^ 0 in Eq. (1281) one can 
immediately see, that the right hand sides of Eqs. (ll5aD - (115bl) depend explicitly on'd, breaking 
rotational symmetry in the x a direction, i.e. 

dt{A a ,C%}± 0. (29) 


It is easy to see that the condition to ensure rotational symmetry of solutions of the field 
equations (TT31) can be written as 

[M, N] = 0 , (30) 

implying that the twisting matrix is diagonal. This happens whenever M does not contain 
terms proportional to a 1 , a 2 . Obviously, the anisotropy generated by twisting matrices not 
commuting with N, would simply rule out the possibility to consider rotationally symmetric 
configurations. Remarkably in the case of a light-like twist vector, due to the decoupling of 
the field equations in the (x 1 ,^ 2 ) plane from the equations with x a = (x°,x 3 ) components, 
(jl5a[) (115b[) . This decoupling allows for solutions which actually break rotational symmetry in 
the x a direction. Whenever M does not commute with N the corresponding M-twisted strings 
have rotationally symmetric spatial sections in any plane orthogonal to the x 3 -axis; however, 
the complete configuration is not rotationally symmetric in the whole space-time. 

Keeping the possibility of breaking rotational symmetry in mind, we now present our Ansatz. 
By singular U(2) gauge transformations (linear in i9) on the scalars in Eq. (1281) we can achieve 
ni = mi = 0. Furthermore, by assuming that the functions, are all real, one reduces 

the number of scalar fields from eight to four (minimality of the Ansatz). Then from Eq. (I22bl) 
it follows that C'J = const, which can be set to zero. 

Finally choosing the radial gauge, our Ansatz can be written as 


A n 


Cl 


A(r,i9)oj a , A r 
C a (r , iH)uj a , C a r 


0 , An = a(r ), 

0, C% = (ci(r),0,c 3 (r)} , $(x‘) 


/</>i (r) 

\(f >2 (r) i> 2 \r)e im J 


(31a) 

(31b) 


We can now display the angular momentum, J for our Ansatz (T3T1) in a more explicit form. 
Exploiting Eqs. (|15aJ) - (ll5bp we obtain 


J = coo (A + Q) , Q = Tr [$ f $(NM + MN) - 2$ t C$N] (32) 

where A = di(adiA) / gf+di(CnDiC a ) / g\. It is worthwhile to point out, that Q is a combination 
of the flavour charge densities [see Eq. (1571) ]. Q = — 4'n a /tg/o;o, with the coefficients n a = 
Tr(iVcr a )/2, n° = Tr 1V/2. 
















Let us now write out the explicit form of the first-order equations (1221) for our Ansatz: 


a' ± yr (0? + 02 + 01 + 02 - 2) = 0 , 

(33a) 

4 ± y' r{<j>\ - 02 + 01 - 0 2 ) = 0 , 

(33b) 

r0i ± ^[(a + c 3 )0i + ci0 2 ] = 0 , 

(33c) 

r02 ± i[(a - c 3 - 2iV)02 + ci0i] = 0 , 

(33d) 

4 ± Q 2 T (0102 + 0102) = 0 , 

(33e) 

r 02 ± 2^ a _ C3 ^ 2 + Cl = 0 , 

(33f) 

t0i ± i[(a + c 3 - 2iV)0i + ci0 2 ] = 0 , 

(33g) 


where for convenience we have chosen units such that £ = 1. The vacuum manifold of S\ for 
the Ansatz (13T1) corresponds to the fix-point manifold of Eqs. (I33af g). a curve, which can be 
parametrized as 

0i = cos a , 02 = sin a , 0i = — sin a , 02 = cos a , 

a = N , ci = — iVsin(2a), C3 = —iVcos(2a:), 0 < a < ‘In . (34) 

3 Twisted vortices 

3.1 Twisted elementary vortices 

The simplest non-Abelian vortex solution is a ‘diagonal’ one, with just 0i, 0 2 , a, C3 being 
non-trivial and subject to Eqs. (j33af g). while 0 2 , 0i, Ci are all zero. In this case the ‘vacuum 
angle’, a = 0. Such solutions have been thoroughly investigated in Ref. [8]. A larger family 
of ‘elementary’ vortex solutions of Eqs. (I33af g). with 02, 01, Ci being also nontrivial, can be 
obtained by a ‘colour-flavour’ transformation from a “diagonal” one [8]. We note that for non¬ 
diagonal NAVs the parameter a is different from zero. For diagonal vortices the relative winding 
between the two doublets is always trivial, N — 0, which remains so for the general elementary 
vortices. The general form of the elementary solution can be written as: 

$ = 0+1 + 0_nV a , A# = a(r ), C# = n a c 3 , 0± = (0i ± 0 2 )/2 , (35) 

where n a is the ‘orientational’ unit vector of an elementary NAV, which can be parametrized 
by the two spherical angles as n a = (sin a cos 0 , sin a sin 0 , cos a). In this case the two moduli 
parameters are just the angles, a and 0 [SUE]. 

We note here that for the most general (with the two moduli) elementary vortex solution 
in Eq. (1351) Cj 0 0, and hence for 0 0 0 it is not in the form of the minimal Ansatz d3T]) : this 
has no influence, however, on our twisting of this solution. It is now simple to generalize, or 
deform the elementary vortex solution by M. Parametrizing the twisting matrix as in Eq. ([7]), 
a short computation shows that the Ansatz 

A = m° , C a — (m-n)[l — C(r)]n a + C(r)m a , m-n = m a n a , (36) 
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containing just a single radial function, C(r), reduces the Gauss constraints, Eqs. (j!5a|) - (jl5b|) . 
to a single second-order inhomogeneous equation for C(r): 

i(rC')' - ic = g\ [0\(C - 1) + i 2 _(C + 1)] . (37) 

Without going into more rigorous mathematics one can easily convince oneself that Eq. (j57j) 
admits a unique solution subject to boundary conditions guaranteeing regularity at r = 0 and 
at r —y oo. Conversely, using maximum-principle-type arguments it is not difficult to show that 
all globally regular solutions are necessarily of the form of the Ansatz in Eq. ()36lh 

The charge density Q, determining the energy and momentum, of the elementary string, is 
given as: 

Q = (m-m - ( n-m ) 2 ) [4>\ + - 2C(r)(j) 1 ^ 2 \ , (38) 

As one can easily see from Eq. (I32j) the angular momentum density of twisted elementary 
vortices vanishes. As already pointed out, this subsection reproduces and extends previously 
obtained results of Ref. |12] using different methods, although our starting was point rather 
different. 

3.2 TWisted coincident composite vortices 

Next we consider a vortex solution with nonzero relative winding between the two flavours. 
Assuming the minimal Ansatz, the vortex in the (a: 1 , x 2 ) plane is rotationally symmetric, sat¬ 
isfying the first-order equations fl33af g) . Such solutions have already been analyzed in Ref. [7], 
where it has also been pointed out that they correspond to superimposed vortices on top of each 
other. Therefore such vortices can be considered as composed of elementary ones. As we have 
already argued, deforming composite vortices with a general matrix, M induces a nontrivial 
angle dependence in the (x° , x 3 ) plane. 

Without losing generality one can parametrize the twisting matrix as 

o ryi 

M = - (l - a 3 ) + — (cos(yu)cr 1 + sin(/i)a 2 ) . (39) 

Taking into account the possible angle dependence of the uj a components of the gauge fields, 
we introduce the following decomposition: 

A = sA 0 + mA + e i{m+ ^ + mA_e~ i{Ni>+ ^ , C a = sG 0 a + mC a + e i{N ^ + mC a _e~ i{Ni}+ ^ , (40) 

together with the conditions A + = A*_, = (C“)* ensuring the reality of the fields; moreover 

Ao, A±, and Cq , C± are functions of r. Then the corresponding Gauss constraints can be put 
in the form 


A r A 0 = gl{rjA 0 + g a C^ - ^) , 

(41a) 

f ab c _ 1 

A -°o + = d(>)CS + V A, - - 4 ), 

(41b) 

A l ''>A= 9 JM ± + D a C|-iV), 

(41c) 

f ab i 

A m C“ + ± 2 iNCl) = gll^Cl + r,‘A ± -- X 5 ), 

(41d) 

A Cn Jcl ± 2i f Cf > - gKgC'l =f j X 2 ). 

(41e) 
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where the SU(2) gauge helds have been split as C a = {C a ,C 2 }, a — 1,3; e 13 = —1, e 31 = 1; 
A r is the radial part of the two-dimensional Laplacian, A (N) = A r — N 2 /r 2 . Furthermore we 
introduce the notations a = { 0 , a} and a 0 = 1 , to present more compactly the often appearing 
combinations 


C (13) = C3C 1 - C1C 3 , 


Va = no sum over A, (42) 

4 = (b?±bf)/2, 


moreover to simplify the formulae somewhat we write rj = 77 ^, r] a = 77 “. We have omitted 
the equation for Cq from Eqs. ill l;ili - il 1 leli . since a straightforward application of the maximum 
principle leads to Cq = 0. The reason behind Cq = 0 is the minimality of the Ansatz (l3Tj) . As 
one can see, the Gauss-type Eqs. il 1 lall - il 1 leli can be decomposed into three equation groups, 
one for {A 0 , Cq}, and one for each {A-t, C±}, decoupled from each other. In fact it is sufficient 
to consider only one of the set of eqs. for ±-components and impose reality on the solutions. 

From Eqs. (j4 1 aft - (j4 left one can easily deduce the asymptotic r —> 00 behaviour of the oj a 
components of the gauge potentials. We note first that for r —y 00 

X° 0 , x 1 ->• cos(2a), x 2 ~i , X 3 ->■ - sin(2a), (43) 

r]\ -A- —r]\ — > sin(2o;), r]\ -A- —rfe — > cos(2a) , 0 , rf A -A- 1, (44) 

then one Ends that for r — y 00 

A) -*• 2 ’ 0 ) Q> 2^2 j C± jX°', C± -A- =F- . (45) 


Quite interestingly, the equations for the angle-dependent components, (14left - (14left can be re¬ 
duced to a quadrature, i.e. to solve a single first-order, linear ordinary differential equation. 
The key observation is that A± = 0 is a solution of Eq. (14left . This is not completely obvious at 
first sight, since assuming A± = 0, Eq. (141cft leads to an algebraic relation/constraint between 
C\. and C±. A straightforward computation shows that this constraint is compatible with the 
remaining two coupled second-order equations (141dft and (14left . As a matter of fact one can 
find yet another simple algebraic relation among the C±. In conclusion Eqs. ill left ill left admit 
a globally regular solution, which can be given as 


A = o, 

rrj 3 Cl' = ~(c 3 r 7 + AT / 3 )<7± + (csy 1 + iNrj^/x 2 )/4 , 


(46a) 

(46b) 


and in terms of the solution of Eq. (I46bft . the remaining functions C±,C± can be found from 
the following algebraic relations: 


7) 1 Ci + 7) 3 C|-jX 0 = 0, (47a) 

x'f'i ± .\ 2 ei + X 3 C"1 = 0 . (47b) 

It follows from Eqs. (14611471) that CA = (7“ and C 2 = —C 2 *, which implies that C± are real 
and C± are imaginary. 
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9i 

a 

QT 

/Qtot 

Vm 

0.4 

0 

0 

6.283 


0.05 

0.0249 

6.222 


0.785398 

3.118 

3.591 

0.77 

0 

0 

6.507 


0.05 

0.0243 

6.261 


0.785398 

2.966 

3.591 

2.33 

0 

0 

6.222 


0.05 

0.0232 

6.211 


0.785398 

2.683 

3.586 


Table 1: Q]^, for #2 = 1 and different values of gi,a. Note that the magnetic energy is 
E 1 = 4vr and 0.785398 « tt/4. 


Using the solution of the Gauss constraints (T471) . the electric energy density, Q, simplifies 


to 


Q = s 2 Q s + m 2 Q m + msQ ms T ($), 
Q, = [1/2(1 - AO - Cg rg\ , 


Qm 

Qms 


1 q \ 0 ^ 

2C i 1 v x „ A , v-V- 


— X I + 0 3 
m J ?v 

-A 0 x°~C*x s + 2Cl 


(48) 


1 s 
TV 2 

?7 3 


V2 + 


s n 
ViX 


2hs J ’ 


where T(i?) = cos (Nfl + /i). In the angular momentum density J [Eq. (132]) ]. Q takes the form 


Q = sNQ s + mNJ m T(d ), 

Jm = ~ A [ 4 ' x ° x2c + _ ■ 


(49) 


The total electric energy, and the longitudinal and the angular momenta are given as 


e 0 j = (Agf + m 2 g^), p = ^ (Agf 1: + m 2 g^), j tot = uqsnq ™, (50) 


where g^ s = f d 2 xQ m , s - 

We have plotted the radial components of twisted coincident vortices for N — — 1. In Figs. 
ID and [21 the planar components are displayed (the background, planar solution). In Figs. [3] 
resp. 01 the Fourier components of the out-of-plane gauge fields are shown. The charge density 
terms Q s and Q m are plotted in Fig. 0 Their integrals over the plane are given in Tabled] 

In Table [H there is a striking difference between the energy of a vortex string with s ^ 0, 
m — 0 and s = 0, m / 0, although they are counterparts in the sense that they have the same 
frequency. Noting that the vortex with m = 0 is rotationally symmetric, while the one with 
m 7^ 0 is not, only the magnitude of the difference is surprising. 

For an explanation of the magnitude of the above-mentioned energy difference, let us apply 
perturbation theory, expanding the solution in powers of a [see Eq. (J34J)], assuming a <C 1. An 
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expansion of the background vortex as 


0i = 0i O) + a 2 0i 2) + • • • , 

02 = a02 1} + • • ■ , 


01 = a0| 1} + ... , 

02 = 02 O) + a 2 02 2) + • • • , 


a = a ^ + a 2 a^ + ... , 
c 3 = 4 0) + a 2 4 2) + ... , (51) 


can be substituted into Eqs. (133al g). yielding a consistent solution. Note, that in the a 0 order, 
the vortex is always gauge equivalent to diagonal one. If a = 7 t/2 , the configuration can also 
be brought to a diagonal form [Tj. The held components A 0 , A±, Cft, C± are expanded as 


4) = 2 + , 

Cl = aC'„ m , 

Cl = ~ + a 2 C 0 3(2) 


and substituting into Eq. 


yields 


C\ = Ci (0) + a 2 Cl (2) 
Cl = Cf 0) + a 2 Cl {2) 


Cl = aC ± 


3(1) 


(52) 


Qs = a 2 


(0{ 1) ) 2 -2C 1 Vi L> - - cd’MD 




( 2 ) 


,3(2)w ,,(0)x2 


Q,n = -2C} (0) <Vf + i t 


(d 0) ) 2 + (v-f) 2 ' 


+ 


+ 


(53) 


The two orders of a between the leading terms of Q s and Q m in Eq. (1531) explain the magnitude 
of the energy difference between the same planar vortex twisted with the same frequency, either 
with a diagonal or with an off-diagonal twisting matrix. 

Finally, we give some arguments for the stability. The conserved charge Q is strongly 
localised, and therefore small perturbations cannot change its value. Planar vortices, being 
absolute minima of the energy in their topological sector, are stable. If there were an instability, 
it would be expected to manifest itself as an energy-reducing deformation along the z axis. In 
the case of the twisted semilocal vortices of Ref. id. with uj a timelike, such deformations 
indeed exist [21], however, they correspond to the same type of instability as those of ANO 
vortices embedded in a two-component extended Abelian Higgs model [22]. In the present 
case, however, such potential instabilities are absent. The spectrum of the perturbation modes 
of the untwisted vortices are gapped, moreover the planar and off-planar perturbation modes 
decouple. Therefore at least for small values of the twist, it cannot change the sign of the 
otherwise non-vanishing positive eigenvalues. 


4 Conclusions 

In this paper, we have constructed charged, stationary rotating non-Abelian vortex strings in 
a I/(2) gauge x SU(2)f[ avor theory. The scalar fields rotate around the string axis, and they have 
a (matrix) phase depending linearly on x a = (t, z) as 

<h(^) = $00) exp f^Mcu a x a 

which is referred to as a twist. We considered here the case uj a u) a = 0, in which, the planar 
equations decouple from those of the t, z components. The energy contribution due to the 
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twist, and the z component of the momentum are both proportional to the Noether charge 
corresponding to the flavor symmetry generated by the matrix M. 

Adding twist to the coincident composite vortices of Ref. |7], leads to some striking phenom¬ 
ena. These vortex strings carry total angular momentum, unless M is purely off-diagonal. If 
M is non-diagonal, the vortex strings are not rotationally symmetric, although, all their planar 
cross sections are. This is explained by the fact, that the nontrivial realizations of rotations 
and ^-translations act on them non-commutatively. The energy of a solution, which breaks 
rotational symmetry, is significantly larger than that of its rotationally symmetric counterpart. 

The analysis of vortex hair for charged rotating asymptotically anti-de Sitter black holes has 
revealed interesting features particularly in the holographic context of the gauge/gravity duality 
|23| . The interplay between the angular momentum of vortices such as those constructed here 
and the black hole angular momentum could give rise to interesting effects. We hope to discuss 
this issue in a forthcoming work. 


Acknowledgements 

This work has been supported by the grants OTKA K101709 and TET 10-1-2011-0071 and 
MINCYT-Argentina HU/10/01. F.A.S. is partially financed by PIP-CONICET, PICT-ANPCyT 
and CICBA. A.L. is partially financed by The Hellenic Ministry of Education: Education and 
Lifelong Learning Affairs, and European Social Fund: NSRF 2007-2013, Aristeia (Excellence) 
II (TS-3647). 

References 

[1] A.A. Abrikosov, “On the Magnetic Properties of Superconductors of the Second Group”, 
Sov. Phys. JETP, 5 (1957) 1174; H.B. Nielsen and P. Olesen, “Vortex-line models for dual 
strings”, Nucl. Phys., B 61 (1973) 45. 

[2] E. B. Bogomol’nyi, “The stability of classical solutions”, Sov. J. Nucl. Phys. 24, 449, 
(1976). 

[3] H. J. de Vega and F. A. Schaposnik, “Classical vortex solution of the Abelian Higgs model”, 
Phys. Rev. D14, 1100, (1976). 

[4] H. J. de Vega and F. A. Schaposnik, “Electrically charged vortices in non-Abelian gauge 
theories with Cherrr Simons term”, Phys. Rev. Lett. 56, 2564 (1986); H. J. de Vega 
and F. A. Schaposnik, “Vortices and electrically charged vortices in non-Abelian gauge 
theories”, Phys. Rev. D34, 3206-3213, (1986). 

[5] A. Hanany and D. Tong, “Vortices, instantons and branes”, JHEP 0307, 037 (2003) 
arXiv:hep-th/0306150 ; 

R. Auzzi, S. Bolognesi, J. Evslin, K. Konishi and A. Yung, “Nonabelian superconduc¬ 
tors: Vortices and confinement in N = 2 SQCD”, Nucl. Phys. B 673, 187 (2003) 
arXiv:hep-th/0307287 ; 

N. Markov, A. Marshakov, A. Yung, “Non-Abelian vortices in J\f = 1* gauge theory”, 
Nucl. Phys. B709, 267 (2005) arXiv:hep-th/0408235 . 

[6] M. Shifman and A. Yung, Supersymmetric solitons, CUP, 2009. 


14 



[7] R. Auzzi, M. Shifman, A. Yung, “Composite non-Abelian flux tubes in J\f = 2 sQCD”, 
Phys. Rev. D73, 105012 (2006) arXiv:hep-th/0511150 . 

[8] M. Eto, Y. Isozumi, M. Nitta, K. Ohashi, and N. Sakai, “Solitons in the Higgs 
phase: the moduli matrix approach”, J. Phys. A: Math. Gen. 39 (2006) R315-R392 
arXiv:hep-th/0602170 ; 

M. Eto, Y. Isozumi, M. Nitta, K. Ohashi, and N. Sakai, “Moduli space of non-Abelian 
vortices”, Phys. Rev. Lett. 96, 161601 (2006) arXiv:hep-th/0511088 ; 

M. Eto, K. Konishi, G. Marmorini, M. Nitta, K. Ohashi, W. Vinci, and N. Yokoi, 
“Non-Abelian vortices of higher winding number”, Phys. Rev. D74, 065021 (2006) 
arXiv:hep-th/0607070 ; 

M. Eto, T. Fujimori, S.B. Gudnason, K. Konishi, M. Nitta, K. Ohashi, and W. Vinci, “Con¬ 
structing non-Abelian vortices with arbitrary gauge groups”, Phys. Lett. B669: 98-101 
(2008) arXiv:0802.1020 [hep-th] . 

[9 1 D. Tong, “Monopoles in the Higgs phase”, Phys. Rev. D69 065003 (2004) 

arXiv:hep-th/0307302 ; A. Hanany and D. Tong, “Vortex strings and four-dimensional 
gauge dynamics”, JHEP 0404, 066 (2004) arXiv:hep-th/0403158 ; R. Auzzi, S. Bolog- 
nesi, J. Evslin, “Monopoles can be confined by 0,1 or 2 vortices”, JHEP 0502 046 (2002) 

arXiv:hep-th/0411074 . 

[10] A. Gorsky, M. Shifman, and A. Yung, “Non-Abelian Meissner effect in Yang-Mills theories 
at weak coupling”, Phys. Rev. D71, 045010 (2005) arXiv:hep-th/0412082 . 

[111 B. Collie, “Dyonic non-Abelian vortices”, J. Phys. A42 (2009) 085404 

arXiv:0809.0394 [hep-th] . 

[12] M. Eto, T. Fujimori, M. Nitta, K. Ohashi, and N. Sakai, “Dynamics of non-Abelian vor¬ 
tices”, Phys. Rev. D84 (2011) 125030 arXiv: 1105.1547 [hep-th] . 

[13] M. Eto, Y. Isozumi, M. Nitta, K. Ohashi, and N. Sakai, “Manifestly supersym¬ 
metric effective Lagrangians on BPS solitons”, Phys. Rev. D73, 125008 (2006) 
arXiv:hep-th/0602289 . 

[14] D. Tong, “The moduli space of BPS domain walls”, Phys. Rev. D66, 025013 (2002) 
arXiv:hep-th/0202012 

[15] A. Hanany and D. Tong, “On monopoles and domain walls”, Commun. Math. Phys. 266 
(2006) 647-663 arXiv:hep-th/0507140 . 

[16] E. Abraham, “Charged semilocal vortices”, Nucl. Phys. B399 (1993) 197-210. 

[17] P. Forgacs, S. Reuillon, and M.S. Volkov, , “Superconducting Vortices in Semilocal Mod¬ 
els”, Phys. Rev. Lett. 96, 041601 (2006) arXiv:hep-th/0507246 ; 

“ Twisted superconducting semilocal strings”, Nucl. Phys. B 751 (2006) 390-418 
arXiv:hep-th/0602175 ; 

Y. Brihaye, L. Honorez, “Twisted semilocal strings in the MSSM”, Int. J. Mod. Phys. A23 
(2008) 581-597 arXiv:hep-th/0701141 . 


15 



[18] Y. Brihaye, Y. Verbin, “Superconducting and spinning non-Abelian flux tubes”, Phys. Rev. 
D77 (2008) 105019 arXiv: 0711.1112 [hep-th] ; J. Garaud and M.S. Volkov, “Super¬ 
conducting non-Abelian vortices in Weinberg-Salam theory - electroweak thunderbolts”, 
Nucl.Phys. B826 (2010) 174-216 arXiv: 0712.3589 [hep-th] . 

[19] L.A. Ferreira, “Exact vortex solutions in an extended Skyrme-Faddeev model”, JHEP 
0905 (2009) 001 arXiv: 0809.4303 [hep-th] . 

[20] P. Forgacs, N.S. Manton, “Space-Time Symmetries in Gauge Theories”, Commun. Math. 
Phys. 72 (1980) 15. 

[21] J. Garaud and M.S. Volkov, “Stability analysis of the twisted superconducting semilo¬ 
cal strings”, Nucl. Phys. B799 (2008) 430-455 arXiv: 0712.3589 [hep-th] ; 
P. Forgacs, A. Lukacs, “Instabilities of twisted strings”, JHEP 0912 (2009) 064 
arXiv: 0908.2621 [hep-th] ; B. Hartmann, P. Peter, “Can type II Semi-local cosmic 
strings form?”, Phys. Rev. D86 (2012) 103516 arXiv: 1204.1270 [hep-th] 

[22] M. Hindmarsh, “Existence and stability of semilocal strings”, Phys. Rev. Lett. 68 (1992) 
1263-1266. 

[23] R. Gregory, P.C. Gustainis, D. Kubiznak, R.B. Mann, and D. Wills, “Vortex hair on AdS 
black holes”, JHEP 1411 (2014) 010 arXiv: 1405.6507 [hep-th] . 


A Notation and conventions 


Here we summarize the definition of the gauge-field strength tensor, covariant derivatives, etc. 
The signature of the flat Minkowskian metric used here is (+, —, —, —), 



d^A,, dpA^, 


dfil - d v C; + e abc Cfc 


c 

V 5 


D^ A =(d, l - l -A„-'-e‘‘cAs> A , (54) 

with a a , (a = 1,2,3) denoting the Pauli matrices, o a o b = 5 ab + ie abc o c . For later use, the 
covariant derivative of adjoint representation fields, D^E 0, = d^E a + e abc C^E c . 

The Yang-Mills-Higgs equations are 


where the color currents are 


d»F, v = g[Jl, 

D»G; U = g 2 2 J a u , 


(55) 


J°= l -Tr[(D^<h-&D,<h] , 

j;= l - Tr [(7^$) V<P - &a a D^] . 


(56) 
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9i 

a 

h 

h 

Pi 

P2 

0.4 

0 

0.07863 

0 

0 

0.6672 


0.05 

0.07855 

0.009484 

-0.01094 

0.6666 


0.785398 

0.05975 

0.1485 

-0.1617 

0.5176 

0.77 

0 

0.1786 

0 

0 

0.9002 


0.05 

0.1785 

0.01635 

-0.01908 

0.8996 


0.785398 

0.1435 

0.2665 

-0.2930 

0.7188 

2.33 

0 

0.4397 

0 

0 

1.2724 


0.05 

0.4397 

0.02963 

-0.03542 

1.2720 


0.785398 

0.4350 

0.5589 

-0.6301 

1.1314 


Table 2: Shooting parameters of planar vortices for r/ 2 = 1 and different values of g\, a 

The flavor current, i.e., the Noether current corresponding to the global SU( 2) flavor sym¬ 
metry is 

= l - Tr , (57) 

where the component has been introduced for the sake of convenience; a U( 1) 

transformation agrees with a gauge transformation with a constant (global) phase. 

B Numerical data 

In this appendix, we present some numerical data of the untwisted vortices (Table [2]) and of 
the twisted strings (Table [3]). The shooting parameters in the Tables are defined at the origin 
as 

</>i (r) = fir 2 + 0(r 4 ) , 

02 (r) = f 2 r + 0(r 3 ), 

and similarly 

A = s 0 o + 0{r 2 ), 

Cl = s w r + 0(r 3 ), 

C'q = S 30 + C(r 2 ), 

where C P = C\ + iC\ and C m = C\- iC\. 

We recall that A± = 0 , Cq = 0, C'f = C’i, C 3 = C 3 , C\ = — C 2 _ and that for a = 0, 
soo = S10 = s 30 = $m = S3 = 0. The relations (1471) give 

-piSm + P 2 S 3 + h = 0 , (60a) 

fiSm ~ /2S3 = 0 . (60b) 

Equations (jUUl) are trivially satisfied for a = 0. For nonzero values of a, the numerical errors 
in Eqs. (1601) vary between 10~ n - 10 -8 . Of the values considered, the minimal error occurs for 
= 2.33, a = 0.05, when Eq. (jfiOb ) is satisfied to a precision of 3 x 10~ 12 , while the maximal 
error occurs for g 2 = 2.33, a = 0.785398, when Eq. flOOh ) is satisfied to a precision of 4 x 10~ 8 . 
Over the intervals shown in the figures, the errors in the algebraic constraints (I47h -bl remain 
below 6 x 10 -5 . 


0 i (r) =pi r + C>(r 3 ), 

Mr)=P 2 + 0(r 2 ) 1 1 j 

C p = s p r 2 + 0(r 4 ), 

C m = s m + 0(r 2 ), (59) 

Cl = s 3 r + 0(r 3 ). 
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,9i 

a 

-Soo 

<SlO 

S30 

Sm 

Sp 


0.4 

0 

0.05 

0.785398 

0.5002 

0.5250 

-0.01637 

-0.2231 

-0.4981 

-0.08282 

-0.0008571 

-0.2007 

0.05893 

0.05881 

0.03250 

-0.007099 

-0.08075 

0.77 

0 

0.05 

0.785398 

0.5005 

0.5682 

-0.02117 

-0.2231 

-0.4980 

-0.08282 

-0.0008310 

-0.1960 

0.09920 

0.09901 

0.05680 

-0.009071 

-0.1055 

2.33 

0 

0.05 

0.785398 

0.5016 

0.7470 

-0.02777 

-0.2926 

-0.4977 

-0.06037 

-0.0007835 

-0.1850 

0.1728 

0.1725 

0.1121 

-0.01163 

-0.1440 


Table 3: Shooting parameters of twisted strings for r/ 2 = 1 and different values of gi, a 



Figure 1: Scalar field profiles of a coincident composite vortex for gi = 0.4, g 2 — 1, a — 0.05 
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Figure 4: Out-of-plane gauge field components, C|, C+, (7+ for r/j = 0.4, g 2 = 1, « = 0.05 
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Figure 5: The electric energy density terms Q s , Q m , for r/i = 0.4, g 2 — 1, a — 0.05 
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